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I. INTRODUCTION 

The electrostatic drift wave possesses the unique feature of being unstable both in the 
fluid domain of a resistive plasma and in the collision-less kinetic description. In both 
cases the electrons play an essential role; in the former case, the instability is due to their 
collisions with heavier species, while in the latter case the instability arises due to the 
Cherenkov type electron interaction with the parallel wave-phase velocity. For plasmas with 
a not so small plasma (/,me /mi ), the mode becomes electromagnetic and in such plasmas 
an inter- play between the electrostatic and electromagnetic part of the mode takes place. 
The electromagnetic component is an Alfvn- type perturbation which is dispersive like the 
kinetic- Alfvn wave (KAW) or the inertial Alfvn wave (lAW). Which of these modes enter 
the perturbations depends essentially on the plasma , yet it also implies a certain range of 
the wave numbers k i 

Shear flows at the plasma edge layers are of critical importance in the development of 
the regimes of high confinement mode (H-mode). The H-mode is defined by a sudden 
suppression of edge turbulence which results in an increase of the edge gradients indicating 
the set up of a transport barrier^. Normally H-mode is accompanied by the appearance of 
a quasi-periodic series of relaxation oscillations involving bursts of MHD activity known as 
edge localized modes (ELMs)^, development of which resulted in a periodic destruction of 
transport barrier and loss of considerable amount of particles and energy from the edge of 
plasma. 

The tokamak edge turbulence is at most weakly collisional. Furthermore, the plasma beta 
is of the order of the electron to ion mass ratio thus introducing electromagnetic effects^. 
Therefore, it is to be expected that weakly collisional drift- Alfven turbulence i.e., drift 
wave turbulence coupled to kinetic Alfven waves, plays an important role in the tokamak 
edge. Obviously, this kind of plasma dynamics should be treated both electromagnetically 
and kinetically^. Despite it is known that poloidal sheared flows suppress the electrostatic 
drift instabilities, the effects of the sheared flows on the small scale electromagnetic insta- 
bilities is still poorly understood. It was obtained in Ref.~, that contrary to the electrostatic 
drift instabilities, the drift- Alfven instabilities in shear flow exhibit the non-modal temporal 
growth. By using the shearing modes (or the so-called non-modal) approach to the lin- 
earised set of the two-fluid MHD equations, obtained in drift approximation for the finite 



beta plasma, it was predicted that in plasma shear flows with hot ions, Tj ^ Tg, the hydro- 
dynamic and resistive drift- Alfven instabilities^ display at time t > {Vq)~ the non-modal 
growth as exp (7oV'o't^) of the amplitudes of the electrostatic potential ip and of the along the 
magnetic field component A^ of the vector potential, where 70 is the growth rate of these 
instabilities in shearless plasma. The developed theory reveals that this growth is temporal 
and ceases at time t > (V^kyPs)' . 



This paper is devoted to the analytical investigations of the processes which are respon- 
sible for the temporal evolution of the kinetic drift-Alfven instability^ in plasma flows with 
E X B velocity shear. This instability is the electromagnetic counterpart of the electro- 
static kinetic drift instability of plasma shear flows, the non-modal kinetic theory for which, 
grounded on the method of shearing modes, was developed in Ref.-. The basic equation 
in Ref.^^ was the integral equation for the perturbed electrostatic potential, which describes 
the multi-time-scales temporal evolution of the instability. The solution^ of this integral 
equation reveals at time t > (Vq)~ the non-modal decrease with time the frequency and 
growth rate of the kinetic drift instability and an ultimate suppression of that instability. In 
the present paper, we extend the shearing modes approach in kinetic theory^ to the kinetic 
theory of the Alfven instabilities of inhomogeneous plasma shear flows. This theory uses the 
transformation of the Vlasov-Amper system for potentials if and A^ to sheared (in spatial 
and velocity) coordinates convected with shear flow and accounts for by this mean the effect 
of spatial time- dependent distortion of plasma disturbances by shear flows. The governing 
equations of this kinetic theory are derived in Sec. II. They compose the system of coupled 
integral equations, in which velocity shear reveals as a time-dependent effect of the finite 
Larmor radius. In Sec. Ill, we applied these equations to the investigation of the tempo- 
ral evolution of the kinetic drift-Alfven instability of plasma shear flow. We obtain that 
contrary to the electrostatic drift kinetic instability, the considered drift-Alfven instability 
experiences the non-modal transient growth with time-dependent growth rate. The renor- 
malized nonlinear theory of that instability, which accounted for the turbulent scattering of 
ions by the ensemble of the sheared drift-Alfven perturbations is considered in Section IV. 
A summary of the work is given in Conclusions, Section V. 



II. BASIC EQUATIONS AND TRANSFORMATIONS 

We use the Vlasov equation for species a {a = i for ions and a = e for electrons), 
immersed in crossed spatially inhomogeneous electric field, Eq (f ) and homogeneous magnetic 
field B II e^, 

^ + v^ + — fEo(r) + Ei(r,t) + l[vx(Bo + Bi(r,t))]')^ = (1) 

ot or rria \ c J ov 

in a slab geometry with the mapping (r, 6, ip) — )■ {x, y, z) where r, 6, ip are the radial, poloidal 
and toroidal directions, respectively, of the toroidal coordinate system. We consider the 
case of plasma shear flow in linearly changing electric field Eq (f) = {dEo/dx) xe^ with 
dE^/dx = const . In that case 

r f)E 

Vo (r) = Vo (x) By = -^-Q^^^y = K^^y (2) 

with spatially homogeneous, Vq = const, velocity shear. The spatially homogeneous part 
of shear flow velocity is eliminated from the problem by a simple Galilean transformation. 
The fluctuating electric field Ei = —Vip — (1/c) dA^/dt, and the fluctuating magnetic field, 
Bi = V X A^ in Eq.([T]) are determined with Poisson equation for the electrostatic potential 

A (f (r, t) = -An ^ Ca fa (v, r, t) dv^, (3) 

a=i,e 

and Amper law for the along magnetic field Bq component A^ of the vector potential A, 

A A^{r,t) = X^ ^c, / f^/a(v,r,t)rfv„, (4) 

a=i,e 

where fa is the fluctuating part of the distribution function Fa {fa = Fa — Foa, where F^a 
is the ensemble average part of Fa). 

It was obtained in Ref.-, that the transition in the Vlasov equation from velocity v and 
coordinates x, y, z to convected coordinates v in velocity space, determined by 

Vx = Vcc, Vy = Vy + VqX, Vz = Vz (S) 

and to sheared by flow coordinates x, y, z in the configurational space, determined by 

X = X, y = y + VQtx, z = z, (6) 



(it is assumed that inhoinogeneous electric field, and respectively shear fiow, originate at 
time t = t(o) = 0) transforms the linearized Vlasov equation for fa = F^ — Foq,, with known 
equilibrium distribution Fo^, to the form, which is free from the spatial inhomogeneities 
originated from shear fiow (see also Eq.(8) in RefJ). In the case of the electromagnetic 
perturbations considered here, that equation is 

+ Vax^r- + [Vay " VaxVot) — — + VazT, ^ UJcaVayT, (^ca + Vq) '"o 
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dA, dA, ^,, dAA dA 



dt "^ \ dx ^ dy ) °^ dy 



ay 

dFoa 



dVn 



VazfdAz ,dAz\dFQa e^ Vgz dA^ dFpg 



nia c \ dx dy J dvax "^q c dy dvay 

{uc is the cyclotron frequency of ion (electron)). Now, with the electrostatic potential (p (r, t) 
and vector potential A^ (r,t), determined by Fourier transforms over the coordinates x, y, 

z as 

ip {x, y,z,t)= [<p {k,, ky, kz, t) e"'^''+''yy+"'^'dk,dkydkz, (8) 



we can obtain the equation for the separate spatial Fourier harmonic of fa- With coordinates 
X, y, z of the laboratory set of references, transformation ([8]) has a form 



^ (x, y,z,t)= ip {k,, ky, kz, t) e''''^^+My-Kt-)+^>'^^dk,dkydkz 



It follows from Eq.([9]), that separate spatial Fourier mode in convected-sheared coordinates 
becomes a sheared mode with time dependent wave number k^ — V^tky in the laboratory 
frame. 

It follows from Eq.([7]), that the transformation of the Vlasov equation to convected- 
sheared coordinates (|5]), (|6]), in which separate spatial Fourier harmonic may be considered, 
converts the spatial inhomogeneity into the time inhomogeneity. That equation reveals 
(see, also, Ref.=2), that the ordinary modal solution to Eq.(I7]), for which ordinary dispersion 
equation may be obtained, exists only at time t ^ (Vq)^ . For larger time, the initial 
value problem, which includes also other time scales, (V^')~ and (Vf^kyPi)' , have to be 
solved. The solution of such initial value problem for the kinetic drift instability reveals. 



that electrostatic potential becomes non-modal^ with time dependent frequency and growth 
rate at time t > (V^y . 

With leading center coordinates X, Y, determined through the shearing coordinates by 
the relations 



x = X 



v± 



■ sm( 



y = Y 



vl 






COS0+1/ot (X-x) 



where rj = 1 + V^jujc-, and with velocity coordinates wj_, 



V^ = V^ COS ( 



^/rjv± sin i 



(10) 



(11) 



determined through the velocities Vx, Vy in convective frame, the Vlasov equation ([T]), in 
which species index is suppressed, transforms into the form 
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(12) 



in which any time dependent coefficients are absent. 

It follows from Eqs. OlOl) . (TTTj) . that in shearing-convecting coordinates a particle gyro- 
motion is different from the gyro-motion determined in convective coordinates flTTl) only 
with spatial laboratory coordinates x, y. Now it consists in the rotation with modified gyro- 
frequency and stretching of gyro-orbit along coordinate y with velocity —Vq {x — X), which 
is negative for x > X and is positive for x < X. With leading center coordinates OlUp the 
Fourier transforms ([8]) for the electrostatic potential becomes 



ip{x,y,z,t)= ip{kx,ky,k^,t)exp 



ik,. X - 
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(13) 



where 



kl (t) = {kx - Vltkyf + -k 



7] 



S/' 



(14) 



and tan 9 = ky/ ^{k^ — V^tky) with the same presentation for A^ (x, y, z, t). It follows from 
Eq. flTSl) that finite Larmor radius effect of the interaction of the perturbation with time 
independent wave numbers k^-, ky, k^ with ion, Larmor orbit of which is observed in sheared 
coordinates as a spiral continuously stretched with time, appears identical analytically to the 
interaction of the perturbation with wave numbers kx — V^tky, ky/y/rj, k^ with ion, which 
rotates on the elliptical orbit that is observed in the laboratory frame. The time dependence 
of the finite Larmor radius effect is the basic linear mechanism of the action of the velocity 
shear on waves and instabilities in plasma shear flow-. 

In what follows, we consider the equilibrium distribution function Fjq as a Maxwellian, 



nojX) 
{2nv^f^ 



exp 



(15) 



with the inhomogeneity of the density of plasma shear flow species on coordinate X. In 
this paper we assume, that velocity shearing rate Vq is much less than the ion cyclotron 
frequency Ud, and put rj = 1. With this approximation we have in Eq. flT2|) that 
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and the solution of the Vlasov equation for the perturbation / (t, X, F, 2;i,t>^, 0, t>^) of the 

distribution function F, f = F — Fq with known Fq becomes 
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(16) 



Using solution ( !T6|) for all plasma species in Poisson equation ([3]) and in Amper law (jlj), we 
obtain the following set of integral equations, which governs the temporal evolution of the 
separate spatial Fourier harmonic of the electrostatic potential if (k, t) and electromagnetic 
potential A^ (k, t) in plasma shear flow, which is capable of handling linear evolution of 
Alfvenic instabilities of plasma shear flows: 
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P (t) (/^ (k, t) = V ^ V / dtJn (h (t) h {t 
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to 
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47r ^ CaSUa (k, t, to) 



(17) 



and 



InEq. (inD 
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I / rftie-5'='"^=(*-*i)' 



to 



'^^^^^ - ky^aeA-. (k, ti)) (1 - klvl^ it - t,f) 
- (l + ikyVde (t - ti) - fc^f^g (t - tif) ck^ip (k, ti 



47r ^ CaSna (k, t, to) = Stt^ X] ^^ / '^^^e' 
A;_L (t) w_L 



■ik^Vzt 
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dv±v±Jo 



UJr. 



fa (t = to, k, W_L, ^^2) , 



(18) 



(19) 



and fa (t = to,k, f^,^^) is the initial, determined at t = to perturbation, assumed here as 
not dependent on 0, of the distribution function Fa. It follows from Eq. (TT9|) . that initial 
perturbation </? (k, t = to) of the self-consistent electrostatic potential is equal to 



V?(k,t = to) 



klito) + kl 



^ eaSnaik,to,to). 



(20) 



In Ref.- we found as more convenient and transparent for further analysis to use the 
exphcitly causal representation of the equation for the electrostatic potential with function 
$ (k, t) = if (k, t)Q {t — to), where 9 (t — tg) is the unit-step Heaviside function (it is equal 
to zero for t < t^ and equal to unity for t > to). The equation for $ (k, t), which is relevant 
for the analysis of the low frequency electromagnetic drift type perturbations, has a form 

t 
{kl (t) + kl) $ (k,t) = --^ [dt,^ {$ (k,ti) [1 - Ao. (t,ti)]} 



to 



t 



+ \2- dtiAoi (t, ti) kyVdi^ (k, ti) + 



to 



t 



where 



+ -^ [ dt,e~'^'''"T.it-tif I ^kyv,, + ,ky^^ (i _ t,)) $ (k, ti) 

■^De J 

to 

■-kzvl, {t - ti) (i — ^^ - ky^deA, (k, ti) U + 47r ^ ejn^ (k, t, to) , (21) 

Ao. (t, ti) = /o (k^ (t) k^ (ti) pf) e-i^K^iW+^i(*0). (22) 

It follows from system ( ITSl) . (I2T1) . that velocity shear meets only in the ions term Aqi (t, ti) 
of Eq.( TT8|) and in time dependent k (t) of the left hand-side of Eq.( l2T1) . that introduce the 
non-modal effects in the temporal evolution of the potentials. 

III. KINETIC DRIFT- ALFVEN INSTABILITY IN SHEAR FLOW 

In this section, we consider the temporal evolution of the kinetic drift- Alfven instability 
in plasma shear flow. Drift-Alfven instability involves the coupling of the density-gradient 
driven drift waves with Alfven waves, when the phase velocities of the two waves become 
comparable. It is obvious, that for this instability it is impossible to obtain the explicit 
analytical solution to system flT7|) . f l?T]) . as well as it was for more simple electrostatic 
kinetic drift instability of plasma shear fiow^. For the approximate solution of the system 
( TT7|) . (EH), we consider the perturbations for which the potentials $ (k, t) and A (k, t) change 
over time as 

\k^VTi\ < |ciln$(k,t)/cit| ~ |ciln A (k, t) /cit| < \k^VTe\ , (23) 



that corresponds to the model of the hydrodynaniic ions and almost adiabatic electrons in 
plasma with /3 = AnnTe/B^ (the ratio of thermal to magnetic pressure) in the range 1 ^ /3 3> 
Trie/mi. Introducing the time scale t* of the order of |(iln$ (k,t) /dt\ and |(iln A (k,t) /dt\, 
the condition 0231) will hold at time ti in Eq.f l^T]) 

1 



^^»it-hr-tt»—^. (24) 

For that time, we have klv'^^ {t — ti) ^1 in electron terms of Eqs.( fT8ll and (I2T1) . so that 
exp (— fc^fj^e (t — ti) ) is not small only for ti ~ t. In ion terms, we have fc^f^j (t — ti) ^1, 
that correspond to small ion Landau damping. Using the approximation e~2'^^^^**^*^*^-* ~ 1 
and using the Taylor series expansion for both potentials in electron terms with accounting 
for the terms of the order of (ti — t) , the system f fTSj) and fl2T]) reduces to the equations 



[-(l + T) + Ao.(t)]$(k,t) 
t 

+ ikyVdi / dti^ (k, ti) A^i [t, ti 



to 



pTT kyVde ^ /, ^ ^ /tT 1 9$ (k, t) 
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2 /C^f Te V 2 A^zf Te dt 



and 
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TT 1 1 52 f.dA{k,t) 
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kyVde^f. Rr ck^ 1 d'^^(k,t) 



'i-m/ttt^::^ ^7^- (26) 



^^"^1 V V 2 fc^f Te / dt 

respectively. Here T = Ti/Te, Vdi{e) = {cTi(^e)/eB) dlnuio/dx is the ion (electron) diamag- 
netic velocity, and 

Aodt) = lo {klit) p^ e-^r^'-^'\ (27) 

In this paper we consider the case of shear flow with small velocity shear Vq, for which 

\kzVA\ > \kyVde\ > IV^)'! ~ 7o, (28) 

10 



where 70 is the modal growth rate of the kinetic drift- Alfven instabihty of the shearless 
plasma. In that case, the non-modal processes, which are determined in Eq. fl25l) by the time 
dependent k± (t) in Aqj (t,ti), are slow and may be determined by the approximation 

—i f uj{]<.,ti)dti — « J aj(k,ti)dti 

$ (k, t) = $ (k) e '0 , A (k, t) = A (k) e *« . (29) 

The second term in Eq.( l25l) . that contains the integral of $ (k, t) over time may be presented 
by the integration in parts in the form 

t H 

/-i J uj{l<.,t2)dt2 
dti^{k)e '0 Aoi{t,ti) 

to 



«l>(k)Ao.(t)e *o +^L^$(k)Ao.(t,t 



a;(k, t) cj(k, t, 



oj 



j'-Oj 



t 



+ ,,„,*(k)yrfi,-j^^_j^e .. j;^(-4«.(M0). (30) 

For exponentially growing potential $ (k, t), i.e. when Imu (k, t) > in (12^ . the first term 
in Eq. (l30!) at time t — to> {Imu (k, t))~ is exponentially greater with respect to the second 
term. Also, the third term in the right-hand side of fl30l) is in \kyVdet\ ^ 1 times less then 
the first one for the times t ^ {kyVde)^ ■ Therefore we retain in the right- hend side of 
Eq. fl30l) only the first term. As a result, system (l25l) . ( 126|1 reduces to the following system 
for amplitudes $ (k) and A (k): 



1 + T - ( 1 - ^\ Ao. (t) + .T,/^-(^ ('^' ') - ^y""'^ 



cu (k,t) J V 2 fczf T. 



$(k) 



- tMM^M^ fl + ^±^i^ A (k) = 0, (31) 

- A;i (t) p^A (k) + CO (k, t) {00 (k, t) - fc,.;,e) ^ (l + ^\/|^^^) ^ (k) 

-||(.(k,t)-V.e)(l + .y|^)^(k)^0, (32) 

which gives the following dispersion equation for the time dependent frequency u (k, t): 



Uji}i,t)) V V 2 k.VTe ) 

_ u;^(k,t) / _ kyVdi \ (1 - ^oi (t)) 
klvl \ u;(k,t)y' ^ 



1-1 / 1 "'j/'^di 



kl it) Pl, 
+ T-Ml-^j^)(l-^o.(t)) = 0. (33) 
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Eq.( l33|) is valid for arbitrary values of A;^ (t) p^^ provided that u (k, t) has a positive imagi- 
nary part, i.e. for the unstable drift- Alfven waves. In long wavelength limit k^ (t) p^^ <^ 1, 
Eg. (1551) reduces to the form 



w 



(k, t) - kyVde) {^"^ (k, t) -u (k, t) kyVdi - klv\) 

TV u (k, t) 



^z'^A^iPs (^ (k, t) - kyVdi) 1 - « 



(34) 



2 kzVTe 

where p^ = p\,^T^^. In the case, when the drift wave frequency kyVde is much less then the 
Alfven frequency kzVA-, Eq.f l34p determines the well known frequency Ui and the growth rate 
7i of the kinetic drift instability, 

1 - kl (t) pI, 



uji (k, t) = kyVd, 



y'-'de 



1 + kl (t) p 



2 ' 



7i(k,t) 



n2;.2 



r^{kyVdeYkl{t)pl{l-kl{t)pl;] 



1 + 



T,; 



2 MTe (l + fci(t)p2)3 V" ■ T, 

and the frequencies a;2,3 (k, t) and damping rate 72 (k, t) of two kinetic Alfven waves 

kl{t)p^ 



(35) 



a;2,3 (k, t) = ±KvA (1 + A;i (t) p^) / , ^^ 3 (k, t) 



TX klv\ 



8 k.VTe {l + kl it) 



2^1/: 



17^- (36) 



It is obvious, that for the damped Alfven waves ( l36l) the approximation ( l30l) is not valid and 
the initial value solution of the system f l25|) . f l26|) is necessary for the proper description of 
the damped perturbations. In the short wavelength limit fcj^ {t) p^ ^ 1, Eq. (!33|) has solution 
for the unstable drift kinetic instability. 



wi(k,t) 
7i(k,t) = 



l^yVdi 



y^de 



27Tk_i (t) pTi 

{kyVdef 
2kzVTekl_ (t) Pi 



T. T,; 



l + 7^ + 



p7;2 



T, nkl{t)k^^c^Xl^\ 






p7;2 



(37) 



and two solutions (valid for kzVAk±pi <^ Wcj) for the frequency and the growth rate for two 
damped Alfven waves, 

'/' [^klv\kl{t)pj 



0)2,3 (k, t) = ±kzVAk± (t) Pi 1 + — 



T. 



72(k,t) = - 



(38) 



Ti J V 8 kzVTe 

So, in both cases of the disparate frequencies of the drift and Alfven waves, Alfven waves 
are damped. In the case when the frequencies of the density-gradient driven drift waves and 
Alfven waves becomes the almost equal, the kinetic drift-Alfven instability develops^. In 
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times t <^ts = (V^kyPi) the right hand side of (1M|) is small. Omitting the right-hand side 
term in (1311) we have solutions 



-i,2(k,t) = ^±(^ + A:M)''', (39) 

which define two Alfven waves, and solution 

^^3 = ky'^de, (40) 

which defines the electron drift wave. In the vicinity of the crossing of the solutions Ui and 
Us the dispersion equation (15B]) becomes 



{6u:{Kt)y{co,-co,) = kiviki{t)pi{co, + \u:,\)i^l-z^-^J, (41) 

where 6u (k, t) is the frequency mismatch, determined as u (k, t) = ui + 6u (k, t) = U2 + 
6uj (k, t). It gives the growth rate of the kinetic drift-Alfven instability for the time t in the 
interval (V"o')"^ < t < t, and for T^ = T^ 

7 (k, t) = ImSoj (k, t) = x —kyVdek± (t) pn^-- (42) 

V J-^ VTe 

For time t > (V^')~ 5 ^^ ~ ^o'^^y ^^"^ therefore 

7(k,t) ^ W— A;yi;deA;j^PTi— Vjj't (43) 

V J-Z fTe 

The solutions for potentials $ (k, t) and A (k, t) will be of the form 

($ (k, t) , A (k, t)) ~ exp ( -iw (k) t-iRe / (Jw (k, ti) rfti 

+ ^J ^kyVdekyPTi — Vot^ j . (44) 

At time t > tg the right-hand part of Eq. (1341) ceases to be small and the frequency mismatch 
becomes of the order of the Alfven frequency and disparity of the drift and Alfven frequencies 
occurs. At these times drift and Alfven modes evolve separately. The drift-Alfven instability 
transforms into the drift kinetic instability, temporal evolution of which was considered in 
Ref.I. 
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IV. RENORMALIZED THEORY OF THE KINETIC DRIFT-ALFVEN 
INSTABILITY IN SHEAR FLOW 

In Ref.- we find, that the effect of the turbulent scattering of the ion Larmor gyration 
angle by the random electric field of the ensemble of the sheared perturbations with random 
phases rapidly suppresses the electrostatic drift turbulence. Here we explore the role of 
that effect on the temporal evolution of the kinetic drift- Alfven instability. It follows from 
Eq. (l32|) . that A ~ {ckz/u (k, t)) $. Then, for the turbulent electric field Ei (r,t) we obtain 
that 

*^~%W»l (45) 

uA k;, ^ ^ 

therefore the turbulent scattering of the ions will be dominated by the electrostatic turbulent 
electric field and the methodology developed in Ref.-^ is completely applicable here. 

[-(l + T) + Ao.(t)]$(k,t) 

+ ikyVd^ / rfti$ (k, ti) Aoi (t, ti) exp -t^ / di y 
to \ tl 



2 kzVxe ' V 2 kzVTe dt 



where the approximation 



((k (t) 5v it) - k (ti) 6r {t,)f) ^ ^ ((k (t) 64> (t) - k (h) 6<p {h)f) 



Cl 



^2t'Jdt^^ (47) 

ti 

was used. In Eq.( l47|) the function C (k, t) is determined iteratively by the equation 
with 

- |2 I - ,2 ( ^ F^ f^TimeV'"^ .\ 
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The integration in parts performed similar with Eq. (l30|) gives 



t 

I 

to 



ikyVdi / dti^{k,ti) Aoi{t,ti)exp -t^ / di 



,C{k,i) 



P 



ti 



kyVd,Aoiit) -^Ju.{Ktl)dt^ kyVdiAoi{t,to) 



u{k,t)-tC{k,t) 



$ (k) e *o 



+ 



a;(k,to) -iC (k,to 



-$(k) 



+ kyVdi^ (k) / dt 



to 



{u;{k,h)-iCik,to)y 



^ I _i/^(k,t2)dt2-t2/dt%ii \ ^ 



X -— e *o 
dti 



— iAodt,h)) 



(50) 



Retaining in the right-hand side of Eq. (!50|) only the first term, which is exponentially growing 
with time, we obtain the following renormalized form of Eq.f l3ip . in which the scattering of 
ions by the ensemble of the shearing drift- Alfven waves is accounted for: 



l + T 



kyVdi 



T 



oj ik,t) - iC ik,t) 

'bJ (k, t) - kyVde) 



Ao. (t) + iT 



n {uj (k, t) - kyVd, 



y'-'de) 



kzVTe 



$(k) 



KzC 



^-"'Itfl-W- 



System (13T|) . ( 15T|) gives the equation for the function u (k, t). With approximation 

(w (k, t) - iC (k, t))-^ ^ uj-^ (k, t) + iC (k, t) tj-2 (k, t) 



(51) 



(52) 



that equation has a form 



fc^f^ 



(w (k, t) - A;j^t;de) (^^ (k, t) - a; (k, t) kyVdi - k^v^) 

„2 /,\ 2 



ki (t) pi {u (k, t) - kyVd^i) - iC (k, t) 



"y 



u{k,t) 



^_^ vrc.(k,i) 



2 fc^fTe 



The solution for the nonlinear frequency mismatch, Su (k, t), follows from f l53|) 



5a; (k, t) = ifczt;^ 



kl{t)pl-^- 



C(k,t) 



1/2 



2 /vT kyVde 



2 V 2 /c^t^Te 



(53) 



(54) 



Wl -h |fcyfdi| 

where wi is determined by Eq.( l39l) . Eq. (l5^ displays the unexpected result: the turbulent 
scattering of the ions does not suppressed the drift- Alfven instability. Instead, that effect ac- 
celerates the non-modal growth of the potentials with time. The quenching of the instability 
occurs due to violation of the frequency matching condition on time t > {V^kyPi)" . 
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CONCLUSIONS 

In this paper, we develop the non-modal theory of the kinetic drift-Alfven instability of 
plasma shear flow. In this theory, the shear flow reveals as the time-dependent effect of 
the finite Larmor radius in the system of integral equations for the electrostatic potential 
and for along-magnetic-field component of the vector potential. This effect is of principal 
importance for turbulence evolution in plasma shear fiows^. It consists in the interaction of 
ions undergoing cyclotron motion with inhomogeneous electric field of sheared modes, which 
due to their distortion by shear flows have time dependent wave number in the laboratory 
frame of reference. In Ref.- we obtained, that this effect reveals in the non-modal decrease 
with time the frequency and the growth rate of the electrostatic kinetic drift instability 
by shear flow and results in the enhanced suppression of this instability. In this paper 
we find, however, that electromagnetic counterpart of that instability - kinetic drift-Alfven 
instability has principally different behaviour in shear flow. We derive, that due to shear flow 
the growth rate of that instability experiences the non-modal transient growth with time 
and, as a result, the anomalously fast development of the drift-Alfven instability occurs. 
The turbulent scattering of ions by the sheared drift-Alfven perturbations accelerates this 
growth. 

At time t > (V^'fcj^Ps)^ , due to the non-modal growth of the mismatch of the drift 
and Alfven waves frequencies, resulted from time dependent finite ion Larmor radius effect, 
this instability transforms into electrostatic kinetic drift instability and separately damped 
Alfven waves. On that stage of the evolution, electrostatic drift instability experiences the 
linear and nonlinear damping and eventual suppression^ conditioned by the velocity shear. 
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Appendix A: THE DERIVATION OF THE FOURIER TRANSFORMED 
Eq.(7) BY USING LABORATORY FRAME SPATIAL COORDINATES. 

Here we display tliat Eq.([7]) may be obtained without the transformation to the sheared 
coordinates (jlj). After the performing the transformation to convective coordinates ([3]) in 
velocity space, Eq.(IT]) becomes 

+ V— r + VqX— :- + UcaVyT. (Wca + Vo) fo 



ay 



dt ' dv ' '> dy ' ''^ 'dva. ' "" "^dv. 

mo, ' dv iTLa c dx dvax "m-a c dy dVay 



lea 
cnir. 



dA^,r>^9A dA dA, 



Sl^n 



(Al) 



dt dy dx 

After the Fourier transforming of Eq.(Al) over the coordinates x,y,z in laboratory frame 
as 

/« (v, k, t) = y U (v, f , t) e-'^'dk (A2) 

we receive the following equation for /„ (t, k, v 



dfa .rff dfa , ■ ff A r , ^/a / , t//\ "^/q 



dk.„ ^ ^ UVnx. 



kv9 ( k, t ) ^;-^ - i—^kr,v^Az ( k, t ) ^^-^ - i—^kyV^A^ ( k, t 



^a V / (TV mo-c v / ar;^; rriaC V / ar;^ 



+ ^^(^-KK^]^. (A3) 

For the deriving from Eq.( 1A3P the equation which couples /„ with potentials (p and A, of 
the separate spatial Fourier mode, as it is in Eq.(I7]), we have to exclude from Eq. flASI) the 
differential operator —Vfiky^, due to which the Fourier mode of fa appears to be coupled 
with all Fourier modes of the electrostatic and vector potentials and depends on the integral 
of (p and Az over wave-number space. The characteristic equation 

dt = -i (A4) 

V'k 

VgHy 

gives the solution k^ + Vf^tky = K^., where K^ as the integral of Eq. (fTO|) is time independent. 
It reveals that fa = fa\Kx-,ky,kz-,t\ = fa [kx + VQtky,ky,kz,t], i.e. the wave number 
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components k^ and ky have to be changed in such a way that kx+V^tky leaves unchanged with 
time. If we use fc^. = K^—VQtky in Eqs. (IA2( ) and flA3p . we obtain for the electrostatic potential 
the presentation ([H]), and we obtain Eq.(I7j) for /„, with time independent K^ = k^, ky = ky, 
kz = kz- Note, that the separate Fourier harmonic in k space may be extracted, if potentials 
Lp and Az also are determined as the shearing modes, i.e. Lp = Lp{kr^ + VQtky,ky,kz,tj = 
(f {kx,ky,kz,t). The obtained results prove, that the solution of the Vlasov equation for 
times t ^ {Vo)~ in the form of the separate Fourier harmonic may be obtained only in 
convected-sheared coordinates. That solution reveals in the laboratory frame as a shearing 
mode (IE]) with time dependent x-component of the wave number. 
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